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Abstract. We revisit the estimation of higher order corrections to the angular power 
spectra of weak gravitational lensing. Extending a previous calculation of Cooray and 
Hu, we find two additional terms to the fourth order in potential perturbations of large- 
scale structure corresponding to corrections associated with the Born approximation 
and the neglect of linc-of-sight coupling of two foreground lenses in the standard 
first order result. These terms alter the convergence (kk), the lensing shear E-modc 
(ee), and their cross-correlation (ree) power spectra on large angular scales, but leave 
the power spectra of the lensing shear B-modc (f3f3) and rotational (ww) component 
unchanged as compared to previous estimates. The new terms complete the calculation 
of corrections to weak lensing angular power spectra associated with both the Born 
approximation and the lens-lens coupling to an order in which the contributions are 
most significant. Taking these features together, we find that these corrections are 
unimportant for any weak lensing survey, including for a full sky survey limited by 
cosmic variance. 
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1. Introduction 

Weak gravitational lensing provides a powerful way to probe the matter distribution 
of the large-scale structure of the Universe and to measure cosmological parameters 
both through growth of structures and geometrical distance projections (see, [H El |Hj 
for recent reviews). As weak lensing experiments improve in precision, it is necessary to 
improve the accuracy to which lensing statistics are computed. Included in the sources 
of theoretical uncertainties are various approximations used in the canonical calculation 
of the weak lensing power spectra. Two such approximations, the Born approximation 
and the exclusion of lens-lens coupling involving two foreground lenses along the line-of- 
sight, were investigated by Cooray and Hu j3j using analytical methods, while numerical 
simulations have also been used to explore these approximations jTUl H2] . Relaxing 
these approximations introduce sub-percent-level corrections to the power spectra of the 
convergence (k) and the shear E-mode (e). Miniscule contributions to the angular power 
spectra of shear B-mode {(3) and image rotational (u) components are also generated, 
though in the standard calculation, these two components are exactly zero. This has led 
to the use of the shear B-mode component monitoring systematics in lensing studies. 

In Cooray and Hu 4J, the estimation of lensing corrections was done perturbatively; 
they computed corrections to the power spectra which were fourth order in the 
gravitational potential of large-scale structure. In this paper, we briefly review these 
corrections and show that there exist two additional terms of the same order. Mainly, 
these terms decrease the power of the convergence and shear e-mode on large angular 
scales. Neither the shear /3-mode nor the rotational component are changed relative to 
the previous estimate in Ref . jl] . Together, the calculation is complete and we can make 
a strong conclusion that the corrections due to the Born approximation and lens-lens 
coupling are negligible for both current and upcoming lensing experiments. Unless one 
begins to probe multipole moments greater than 10 4 , it is likely that these corrections 
can be ignored for any experiment, including an ultimate all-sky experiment limited 
by cosmic variance alone. While the calculation related to higher order corrections is 
complete, there may be other sources of error in the theoretical estimates of lensing 
statistics [3 El- For example, theoretical estimates of the underlying non-linear dark 
matter clustering is not known to better than 10% at arcminute scales and is likely 
to be further complicated by baryon physics [TTH ITIj . Another source of important 
error is that lensing statistics are based on reduced shear rather than shear directly, 
which results in perturbative corrections from convergence to shear [T5] . Clustering 
of background sources also contaminates weak lensing statistics [Zj and may be an 
important source of uncertainty for tomographic lensing measurements |Hj in which the 
background source redshift distribution is binned to obtain lensing measurements as a 
function of redshift. Based on our calculation, we are now confident that that these 
remaining issues dominate the theoretical calculation and interpretation of the weak 
lensing angular power spectra. 

This brief paper is organized as following: In the next Section, we outline the 
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calculation related to higher order corrections to weak lensing angular power spectrum. 
We identify two new terms by extending the calculation of Ref . [3] , and make numerical 
estimates of the correction. We find that the overall correction is negligible even for an 
all-sky lensing experiment limited by cosmic variance and probing fluctuations out to a 
multipole of 10 4 in the angular power spectrum. We conclude with a brief summary of 
our results in Section |31 For illustrative purposes, we assume a flat ACDM cosmological 
model with parameters Q c = 0.3, fi& = 0.05, Q\ = 0.65 for CDM, baryon, and 
cosmological constant densities relative to the closure density, a scaled Hubble constant 
of h = 0.7, and a primordial power spectrum of fluctuations normalized to erg = 0.85 
with a tilt of n = 1. 



2. Calculational Method 



Here we provide a brief summary of the analytical calculation. For further details, 
we refer the reader to Cooray & Hu jlj. For simplicity, we also use their notation 
throughout. 

The distortion tensor for a weakly lensed source is given by an integral over a 
photon's path, x, through the large-scale gravitational potential: 

rxs 

ip ab (n,x s ) = 2 / dx g(x,Xs)®,ac(x;x)[$cb + ipcb(n,x)}, (1) 
Jo 

where n is the source's sky position, Xs is its comoving distance from us, $ is the 
potential and 5 is the Kronecker delta. Commas denote spatial derivatives in the 
transverse directions, and there is an implicit sum over repeated indices. The weighting 
function g assumes that the source is at a fixed distance away; it is given by 

tU,x) = {' A(X ~* yUW), ' U< X) (2) 
10 for x > X 

with dA as the angular diameter distance. The power spectra of weak lensing observables 
are found by combining components of C abc d, defined as 

(PMMW) = (2*W - iOc«w(i), (3) 

with ^(1) the 2D Fourier transform of ^(n). 

The canonical calculation of the power spectra takes advantage of the fact that the 
potential - and therefore the deflection of light - is small. Lens-lens coupling, i.e. the 
appearance of ip on the right side of equation Eq. [H can clearly be ignored to first order 
in $. Furthermore, since the photon path is nearly a straight line, we may write 

x(n, x) = nd A (x) + <5x(n, x), ( 4 ) 

and in the Born approximation set the transverse deflection Sx to zero. The result is a 
simplified integral over the undeflected path, 

ip ab (n,xs) = 2/ dxg(x,Xs)$A nd A(x);x) + 0($ 2 ), (5) 

Jo 
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which, in the Limber approximation leads to the standard expression for the two 
point functions of ip in terms of the power spectrum of <3? ^7] 

C abcd (l) = UJbUd J d X 9{ *^ 2 P(l/d A ; X ) ■ (6) 
The first order spectra of the lensing observables («, e, /3, to) are 

cr = cr = c? = | d X 9{ ^ 2 p{i/d A - x) (?) 

with all other combinations equal to zero. Note that throughout this paper, we will be 
employing Limber approximation to calculate angular power spectra. This assumption 
makes use of the flat-sky approximation and could potentially affect the large angular 
scale correlations. However, as a test case, we considered the calculation exactly 
by integrating over spherical Bessel functions that one encounters with the all-sky 
projections and found our results to be accurate at the percent level and below at £ > 10. 
Just as the Limber approximation is adequate to calculate convergence power spectrum 
at the first order (e.g., Ref. [TT] based on simulations), we believe it is adequate for the 
calculations at the higher order. 

2.1. Perturbative Corrections 



To relax the Born approximation, we simply Taylor expand the potential in equation Eq. 
^ about the undeflected photon path. With the transverse deflection given by 

dAx) 

d A { X ' 



5x a (n, x) = -2 / d X 'g(x', X)^$f ,a(x; x!) , (8) 
the potential is 



$(x; x) = ®(nd A ( X ) + 5x; X ) = $(nd A ; x) + Sx a $, a (nd A ; X ) 

+ -5x a 5x b $ jab (nd A ; X ) + 0($ 4 ) . (9) 

Note that 5x depends on the deflected path, and must also be expanded perturbatively. 
Substituting equation Eq. El into equation Eq. |HJ we find the deflection to second order 
in potential fluctuations to be: 

5x a (n, X ) =5xV(n,x)+5x£\n,x)+0(<S> 3 ) (10) 

5xM(n,x) =-2 J dx'g(x',x)j^^A^d A (x');x') (ii) 

5x£\n, X ) =-2 J ^(x / ,x)^^$,a fe (n^(x / );X , )^ 1) (n ) x')- (12) 

Now we may express $(x), the potential at the deflected position, using integrals over 
the undeflected path: 

$(x; x) = Hnd A ; X ) + {Sx® + 5x^)$ )a (nd A - X ) 

+ Isx^SxP* , ab {nd A] X ) + 0($ 4 ) . (13) 
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Inserting equation Eq. US] into equation Eq. [T] leads to ip as a line-of-sight projection of 
a source field, 



<Mn,xs) = 2y dx g(x, Xs)S a b(nd A ; x), (14) 
where Sab(nd A ; x) includes the following terms (all depend on x): 

S$ =<S>And A ) (15) 

S% B) =fe«(n)$, af)C (nrf A ) (16) 

S%* ) = 6x® (n)<& a6c (nd A ) (17) 

fig^ = (1/2)*^)**?^) W^a) . (18) 

These terms yield the standard second order expression for C a & c d plus higher order 
corrections. Third order corrections vanish in the Limber approximation. Cooray and 
Hu [3] computed all fourth order corrections except for a term that couples S( 3X > to 
We find that it is given by 

^L(l) = 4 J d ^-^f P ^x) (19) 
rf2/ ' /, 7 i/,/ , , wi i/N2 f i,9(x',x) 2 d A (x) 2 



PLiOi x) = - 4/ 7^r- 2 (W d + I'J'Mii ■ i') 2 / 



xi fe*Mdy)4 (20) 

where -Pj^ is the correction to the power spectrum of the source field. The corrections 
to the power and cross spectra of the lensing observables (k, e, (3, u) are 

x /*wKdrt i *Md*>4 (21) 



(22) 



with 

2 for a/3 = 
3/2 for a/? = ne 
1 for a/3 = ee 
otherwise 

Similarly, lens-lens coupling can be accounted for by iteratively expanding the 
recursive term in equation Eq. ^ This yields the following second and third order 
source terms: 

= - 2^ >ac (nd A ; x) J d X 'g(x', x)®A n Mx'); x') (23) 

S S, L) = 4$,ac(nrf A ;x) / dx'g(x\x)^,cd( nd A(x');x) 

x fdx"g(x",x'^M^d A (x"y,x"). (24) 



The first-second order correction to C a b c d vanishes under the Limber approximation. 
Cooray and Hu |lj calculated the second-second order correction, but neglected the 
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correction involving the first and the third order terms. This term does not vanish; it 
can be written as 

^22 = MX)- 6 / d X 'J d 2 ne^(sf\nd A ( X );x)sS l L \0;x')) 

+ (ab «-> cd) (25) 

= 4d A ( X )- e J jj d X 'd X "d X '"J d 2 ne* n g( X ",x'Mx",x") 

x ($ >a6 (nd A (x); x)$, ce (0; x')$ e/(0; x")$ *"')) 

+ (ab «-> cd) . (26) 

To complete the derivation, each of the four $ s above must be expanded into its Fourier 
modes, 0(k). Then to calculate the ensemble average, (0(k)0(k / )0(k")0(k" / )), one must 
sum over all possible pairs of Wick contractions of the s. Contracting two s ultimately 
produces delta functions which match up the corresponding %s. For example, these 
contractions ^ ^ 

(0(k)0(k')0(k")0(k"')> 

lead to a term containing 5(x — x")^{x' ~ x'")- m fact, this is the only non-trivial 
term: because of the weighting functions in Eq. EEl delta functions in other terms 
{$(x" ~ x') an d S(x" — x'")) wn l cause the integrals to vanish jlj. Nonetheless, the 

(13) 

first-third correction does not vanish completely, and we find that Patcd * s coincident ally 
equal to P^ cd in equation Eq. I2TI1 

There also exist corrections due to couplings between Born and lens-lens source 
terms. However, in this case, we do not find additional corrections apart from those 
already considered by Cooray and Hu [3]. To the highest order in importance, i.e., 
fourth order in potential perturbations, this completes the calculation of corrections to 
weak lensing angular power spectra. 

3. Results & Summary 

We have calculated the 0($ 4 ) corrections to weak lensing power spectra that arise from 
the Born approximation and lens-lens coupling; they are illustrated in Fig. Fig. ^ The 
primary effect of the new term, , is to make the total correction to the kk, ne and 
ee spectra negative for / -C 200. The (3 and u> spectra are unaffected by the new term 
when comparing to the calculation of Ref. |3j. As discussed in the u> term arises 
from lens-lens coupling only; it is not a test of the Born approximation as was suggested 
in Ref. [TO]. The larger difference in the amplitude of the rotational power spectrum 
shown here and the one shown in Ref. |H] (their Figure 2) is due to the large difference 
in the source redshift used between these two calculations: here we use z s = 1 for 
galaxy lensing surveys while z s = 1100 is used in [S] for lensing of cosmic microwave 
background anisotropies. The large difference suggests that non-linear couplings that 
lead to these second order corrections grow significantly with increasing path length 
between source and observer, though the correction still remains below the first order 
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term. The conclusions of Cooray and Hu are unchanged: corrections related to the Born 
approximation and lens-lens coupling can be safely neglected for current surveys. They 
do not limit the use of /3-modes in monitoring systematic errors. 

Furthermore, since our calculation is now complete, we can make even a stronger 
statement. As shown in Figure Fig. ^ the higher order corrections are below the cosmic 
variance level of an all-sky survey, y^2/(2Z + l)Cj, out to a multipole of 10 4 when these 
corrections start to be higher than the cosmic variance level. This suggests that for all 
experiments that use statistical information out to a multipole of 10 4 , the angular power 
spectra will not be affected by the corrections associated with the Born approximation 
and the coupling of two lenses. This, however, does not imply that theoretical estimates 
of the angular power spectra are accurate to such a high precision. Throughout this 
paper, we have relied on the Limber approximation, which ignores contributions to 
the power spectra from modes with wave- vectors parallel to the line-of-sight. This 
approximation is justified by lensing simulations JTj, which agree with the canonical 
calculation of the power spectra and hence our conclusions at the few percent level. 
We have also assumed that the weak lensing calculation is valid, despite the fact that 
the lensing observables, such as convergence, have non-gaussian distributions with tails 
that can affect the non-linear part of their spectra, depending on how one includes high 
lensing peaks that merge weak lensing scales to strong lensing in the data analysis |12j . 
Beyond issues related to the weak lensing approximation, the non-linear matter power 
spectrum for a given cosmology is uncertain at the 10% level, both due to uncertainties 
in the mapping between linear to non-linear dark matter fluctuations and the effect 
of baryons in determining how dark matter clusters on arcminute scales and below 
|13| Hlj . Accounting for such uncertainties are beyond those possible with analytical 
calculations. In fact, all remaining uncertainties are those that must be addressed with 
numerical simulations. We can now be confident that the basic analytical formulation 
of lensing power spectra is well explored and the first order result is accurate enough 
for all upcoming experiments. 

To summarize, here we have revisited the estimation of higher order corrections 
to the angular power spectra of weak gravitational lensing and have extended the 
calculation of Cooray and Hu jlj that first described the corrections related to the 
Born approximation and the neglect of line-of-sight coupling of two foreground lenses 
in the standard first order result. We found two additional terms to the fourth order in 
potential perturbations of the large-scale structure and these terms, in return, altered 
the convergence (kk), the lensing shear E-mode (ee), and their cross-correlation (ne) 
power spectra on large angular scales. There is no modification to the lensing shear 
B-mode ((3(3) and rotational (ujuj) component power spectra when compared to the 
previous estimate. With these new terms, the calculation of corrections to weak lensing 
power spectrum associated with both the Born approximation and the lens-lens coupling 
is complete. The overall numerical result is that these corrections are unimportant for 
any weak lensing survey including for a full sky survey limited by the cosmic variance. 
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Figure 1. Left: Weak lensing power spectrum and corrections. The solid line 
shows the first order spectrum; the dot-dashed line shows the new correction, 
Cf " (77 = 1), which is negative; the dotted and dashed lines show the full 
corrections (Born and lens-lens) to the kk and ee spectra respectively - they 
are negative at low I and positive at high / as indicated by (-)/(+) signs. 
For reference, we also show the cosmic variance level of an all-sky experiment 
with a thin long-dashed line. The comparison shows that for I < 10 4 , the full 
corrections from Born and lens-lens coupling are below the cosmic variance level 
and is unlikely to be an error of statistical significance. Right: Full corrections 
divided by the first order result. The dotted and dashed lines show the relative 
corrections for kk and ee respectively, while the thin long-dashed line is the 
ratio of cosmic variance to the first order result. 
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